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Using renormalization-group methods, we derive differential equations for the all-orders summa- 
tion of logarithmic corrections to the QCD series for R{s) = a{e^e~ hadrons)/(T(e'''e~ fi'^ fj,~), 
as obtained from the imaginary part of the purely-perturbative vector-current correlation function. 
We present explicit solutions for the summation of leading and up to three subsequent subleading 
orders of logarithms. The summations accessible from the four-loop vector-correlator not only lead 
to a substantial reduction in sensitivity to the renormalization scale, but necessarily impose a com- 
mon infrared bound on perturbative approximations to R{s), regardless of the infrared behaviour 
of the true QCD couplant. 



For center-of-mass squared-energy s, QCD corrections to R{s) = cr(e+e hadrons)/o-(e+e fi^^i ) are scaled 
by a perturbative QCD series (S): 
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This series is extracted from the imaginary part of the MS vector-current correlation function 0, |3| , 

oo n — 1 
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with coefficients T„ ,„ tabulated in Table ||] for 3-5 active flavors, as appropriate for the choice of the center-of-mass 
squared-energy s. Each order of this series depends upon the MS renormahzation scale parameter fi, both through 
the couplant 



and through powers of the logarithm 



x{p) = as{ii)/n 
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Nevertheless, the all-orders series S must ultimately be independent of renormalization scale. R{s) is a measurable 
physical quantity necessarily independent of ^, the artificial scale entering QCD calculations as a by-product of the 
regulation of Feynman-diagrammatic infinities. Hence, 



^ dS[x{fi),L{fi)] 



The above renormalization group equation (RGE) is simply a chain-rule relation in which 



p{x) = 
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where known Q MS /3-function coefficients Pk are also tabulated in TableQ] Thus, the RGE ((Sj is generally employed 
to provide scale dependence to the couplant x. 

"Optimal" renormalization-group (RG) improvement is the inclusion of every term in a perturbative series of the 
form ^ that can be extracted by RG- methods from a perturbative computation to a given order For example, 
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TABLE I: CoefRcients for the imaginary part of the four-loop-order vector-current correlation function, as well as coefficients 
for the four-loop order MS /3-function, are listed for three, four and five quark flavors. 



a next-to-next-to-leading (NNL) order perturbative calculation determines only the coefficients Ti g, o? s-nd T2_i of 
the series However, the RGE can be utilized to determine all coefficients Tk+i^k and Tk+2.k within the series 
(0) . The contributions of this infinite set of coefficients may then be summed analytically, as described below, thereby 
providing an "optimal RG-improvement" of the NNL expression. 

The series S[x,L], as defined in Eq. (j^J, may be rearranged in the following form: 

00 

5[x,L] = l + ^x"5„(xL), (7) 

where 

= ^ T„,„_feu"-^ (8) 

n—k 

Given knowledge of the A:*''-order series coefficient T^, g = Sk{Q), one can obtain Sk{xifJ.)L{ii)) explicitly, thereby 
summing over the entire set of fc*''-order subleading logarithms contributing to the series Q. If we substitute the 
/3-function series JHl into the RGE (jSJ, we find that the aggregate coefficient of x"-L"^p vanishes {n > p) provided 
the following recursion relation is upheld: 

= (n - p -t- l)T„,„_p+i - ^(n - £ - l)l3^Tn-e-i,n-p- (9) 

e=o 

For example, if p = 2, this recursion relation [T„.„_i — PoTn-i.n-2\ relates all leading-logarithm coefficients T„^„_i 
within the series to the known coefficient Ti.o — 1, thereby enabling one to sum all orders of the leading-logarithm 
contributions 



,SiixL) = X Vr„,„_i(xL)"-i = (10) 

1 — pnxL 



PqxL 
to the series S[x, L\. 

More generally, the recursion relation ^ may be utilized to obtain a succession of first-order inhomogeneous linear 
differential equations for the functions Sk(u) within Eq. ((JJ). If one multiplies Eq. jnj by and then sums over n 
from n ~ p to 00, one finds from the definition ^ of Sk{u) that 

which can be trivially rearranged (k = p — 1) into a set of first-order linear differential equations 

du (1 - Pou) (1 - Pou) j-^ \ du J 
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with initial conditions Sk{0) = T^^. 

Noting that S'o(u) — 1 and that Ti q ~ 1 regardless of the number of active flavors, we find the first four solutions 
of lO to be 
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To explore the near-infrared regime of perturbative QCD, we specialize to the case of three active fiavors. Using 
TableQln/ = 3 values within Eqs. p3|l - l|l()() . we find that the version of the series {Tjl which incorporates full summation 
of leading and two subsequent subleading orders of logarithms is given by 



S'-^^ [x{p),L{p)] = 1 + xSi{xL) + x^ S2ixL) + x^ SaixL), 



where {u — xL) 
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Si{u) = 1/(1 -9w/4), 
1.63982- if log(l-9u/4) 



1.31057 
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3.16049 log^ (1 - x)} 



{i-luY (1 

Moreover, we note from Eq. (|16|l that 

5.35589 [-6.62811 + 4.65981 log (l - f xL)] 
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[r4,o + 11.9840 + 31.8738 log (l - fxL) + 29.5946 log^ (l - fxL) - 5.61866 log^ (l - fxL)] 
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thereby providing for inclusion of the x'^SAixL) contribution to the series l(7|) for three active flavors. The series 
coefficient T/^^ appearing in Eq. (|21|l has not yet been calculated perturbatively, which is why we have not included 
the x'^Si{xL) contribution to S\x,L] in Eq. p7|l . (An asymptotic Pade approximant estimate T4^o — 1-90 for the 



n/ = 3 case is presented in ref. 

To examine whether the summation of leading and subsequent subleading logarithm factors decreases dependence 
on the unphysical renormalization-scale parameter /x, we compare the /i-dependence of Eq. (|17() for a flxed value of s 
to that of the n/ = 3 version of the series Q truncated after four- loop-order {M) contributions to the vector-current 
correlation function: 



S^^'^^ [x{^i),L{^i)] = 1 -I- a; + (1.63982 + 9L/A)x^ + (-10.2839 + 11.3792L + 81LVl6)a;3 



(22) 
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Such ^-dependence enters Eqs. (|17ll and (|22|l both through L = log(/x^/s) and through x = x{^), which is assumed 
to evolve via Eq. © (with = 3 choices for /Sq-s) from an initial value choice x{mT) — as(m7.)/7r = O.SS/tt j51,l3- 
Figure m displays a comparison of the /i-dependence of Eqs. (|17|l and (|22(l at fixed s = 1.5 GeV^ . Although both 
expressions exhibit little variation with fi over the 1.3 GeV < /i < 3 GeV range, we see that Eq. H17|l exhibits much less 
variation with in the near-infrared regime below 1.3 GeV. These results clearly indicate that renormalization-scale- 
invariance is more effectively upheld via the summations of leading and subsequent subleading orders of logarithms 
that occur within Eq. (|17|l . We emphasize that Eqs. p7l) and l|22|) both follow from "RG- improvement" of the same 
calculational information [the coefRcients Ti g, 72.0; Ts.o]; simply put, such RG-improvement is more effectively 
implemented in Eq. H17|) than in Eq. 
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FIG. 1: The renormalization-scale (fi) dependence of the summation-of-logarithms (E) series (I17II is compared to that of the 
series 11221 , which is truncated after four- loop (4L) contributions to the vector-current correlation function. For both series the 
physical momentum scale s is fixed at 1.5 GeV^. The evolution of the couplant x{jj.) for both series is referenced to the initial 
value as{mr) — 0.33. 

The usual prescription for obtaining the purely-perturbative (non-power-law) QCD contributions to R{s) at four- 
loop order 0,0 is to set ^ — y/s within Eq. 122|l . and then to substitute the resulting series, 

S'('*^) [a;(V^), L(V^)] = 1 + a;(V^) + 1.63982 x^{y^) ~ 10.2839 x^{y/^), (23) 

into Eq. [Note from Eq. Q that L{y/s) = 0.] This prescription follows from the presumed renormalization-scale 
invariance of the truncated series H22|) , thereby leading to an expression that depends only on the physical scale s for 
the electron-positron annihilation process. Note that all s-dependence of Eq. (|17|l resides entirely in the variable L. 
Let us first fix /i = itIt within Eq. H17|l so as to incorporate the benchmark couplant value xirnT-) — O.SS/tt everywhere 
X appears in Eqs. (|18|I - H20|) . With this choice, the following summation-of-logarithms series can be substituted into 
Eq. iP) to obtain R{s): 

o(S) r / ^ Tf M 1 0-33 „ /0.33 , (ml 
S^^> [x{mr), L{mr)] =1 + log ^ 

TT \ TT \ S 

In Figure |21 we compare the s-dependence of this series to that of Eq. (|23|) , for which all s dependence resides in the 
evolution of x(^/s). To make this comparison, such evolution is anchored to the initial value x^mr) = O.SS/tt via the 
differential equation © with Uf — 3 values for /3o-3 [Table This initial value ensures that the series (|23|l and H24(l 
coincide when y/s — rrir- Figure |21 shows that both series continue to coincide over the range 750 MeV < y/s < rrir- 
For values of \/s less than 750 MeV, however, the truncated series (|23|l drops off quite suddenly at \/s ~ 650 MeV, a 
consequence of the large negative coefficient of x^ (a/s), whereas the series (|24|l [as obtained from the full summation- 
of-logarithms series (I17ll ] continues to probe the infrared domain of R{s) even for values of -^s ~ 400 MeV. In short, 
the summation of all leading and subsequent two subleading logarithms within the perturbative series l(2Jl serves to 
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extend the domain of the R{s) series further into the infrared. This property, as well as the reduced renormalization- 
scale dependence evident in Figure ^ suggests that such summation is particularly appropriate for the near- infrared 
region characterizing sum-rule applications of purely-perturbative QCD corrections to current-correlation functions. 
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FIG. 2: The center-of-mass squared-energy (s) dependence of the four-loop (4L-) truncated series 12311 is compared to that of 
the summation-of-logarithms (E) series l|24|l . as described in the text. In both series, as(mr) is taken to be 0.33 so that the 
series equilibrate at ^/s = m^. 



It is also evident from Figure 12 that the domain of the summation-of-logarithms series H24|l manifests a singularity 
below = 400 MeV, despite the fact that the s-dependence of the series (|24f) is decoupled entirely from any infrared 
behavior of the couplant x, which is held constant at a;(mT-). To understand this restriction on the domain of R{s), 
we first note that each summation (|13|I - H16() becomes singular when 1 — (3qxL —> 0. Such resummation singularities 
have also been observed to occur in completely different contexts, including deep inelastic structure functions 0]. 

The singularity property of eqs. H13(l - (|16|l is upheld for all summations Sk{xL). The solution to the differential 
equation H12|) is necessarily of the form 



Tk,o 



Sk(xL) — — ^ — — t: -t- (Particular solution depending on{S'A;-i, Sk-2, <S'i}) • 

(1 - PoxLf 



(25) 



Since the coefficients Tk^ are results of k -order Feynman diagram calculations, the k -order pole in (|25() at 1 — 
(3(jxL = is genuine and will not be canceled by particular-solution contributions that are sensitive to at most 
{k — l)*''-order Feynman-diagrammatic coefficients {Ti g, 72,01 • • • j ^fc-i.o}- For ^ given choice of renormalization 
scale n, this singularity implies [via Eqs. @ and that each summation Sk{xL) within the full series Q becomes 
singular for a sufficiently small value of s: 



1-/^0^ log 



= 



= /i exp 



(26) 



For example, if the renormalization scale /.t is chosen (as in Figure|2Jl to be mr, a choice for which as(mT-)(= 0.33±0.02 
0) is phenomenologically accessible, then each term in the series ((TJ is seen to become progressively more singular 
as s approaches exp[— 47r/ (9 x 0.33)] = (215 MeV)^ from above. Furthermore, Figure O shows that the low-energy 
behaviour of the resummed expression H17|) is only weakly dependent on fj, in the region ^/s > 600 MeV for choices of 
/i between 0.6mT- and l.Gnir, with the appearance of singular points corresponding to 195 MeV < Smin < 250 MeV 
from 

It is to be emphasized that this infrared boundary on the physical scale s entering term-by-term within the series 
Q is not a manifestation of any infrared boundary 9] on the evolution of the QCD couplant x{fi). Even if the higher 
order contributions to the /3-function ((HJ were somehow to conspire to allow the couplant to be well-behaved in the 
infrared region [e.g., to have infrared-stable fixed point behavior], the restriction (|26|1 would still apply upon making a 
specific choice of the renormalization-scale parameter /i and its corresponding value of as(/i). Curiously, though, this 
infrared restriction on s can be easily shown to coincide with the "infrared-slavery" Landau singularity A associated 
with naive evolution of the QCD couplant via a one- loop /3-function. The one loop version of Eq. (0) , 



2 '^^ a 2 '^s 
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FIG. 3: The s dependence of the four-loop truncated series 11711 for equally-spaced values of the renormalization scale in 
the range O-Gm^ < fJ. < 1.6m-r. The couplant x{^) is obtained from three-flavour four-loop evolution via Q from the initial 
condition x (mr) — O.SS/tt using TableQlvalues for the /3 function coefficients Po-Ps- 



is satisfied by the relation 



/3olog(AiVA2) 



(28) 



which is equivalent to Eq. 1)261) provided Smin is identified with A^. Indeed, in a one- loop world [x^^{ii) = 
l//3o log (/i^/A^)] where A serves as a universal infrared boundary, the one-loop analogues of the summation-of- 
logarithms series ((17)1 and the truncated series are necessarily equivalent: 



5(p = l + x(")(M)5l (m) log 



= 1 



/?olog(^) 



l-/9o 



00 log 



log 



(29) 
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We find it remarkable that A, the one-loop couplant's Landau pole, persists as an infrared boundary on the domain 
of each summation contributing to Eq. |(7J), the summation-of-logarithms formulation of the perturbative series within 
i?(s). Consequently, s = A^ serves as an infrared boundary for any approximation to R{s) involving the truncation 
of the series such as the expression H17|) obtained via optimal RG-improvement of the four-loop vector-current 
correlation function's imaginary part. 
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